Does deuteron-induced fission on actinide nuclei prevails over the breakup at low 

incident energies? 
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An analysis of the ^''^Pa(d, Sn)^"^"!] reaction excitation function at energies around the Coulomb 
barrier has taken into account the pre-equilibrium and compound-nucleus cross sections corrected 
for the deuteron-breakup decrease of the total reaction cross section, as well as the inelastic breakup 
enhancement. The analysis reveals the dominance of the deuteron breakup mechanism unlike a 
former assessment in this respect of the deuteron-induced fission process. 
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I. INTRODUCTION 



shown to be valid in the presence of the deuteron breakup 
without assuming specific breakup reaction mechanisms. 
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The analysis of the deuteron-induced reactions at low 
energies, in terms of the usual nuclear reaction models, 
is challenging due to deuteron breakup (BU) following 
its weak binding energy, Bd='2..22A MeV. The various re- 
actions initiated by the breakup neutrons and protons 
render the study and prediction of the deuteron reaction 
cross sections more complex. This is why recent mea- 
surements of the 23ipa(d,3n)230U and 23ipa(p, 2n)230u 
reactions cross sections, between 11.2 and 19.9 MeV T], 
and respectively 10.6 and 23.8 MeV [21], are particularly 
useful for the analysis of breakup effects on the former 
excitation function. The outgoing energy of the breakup- 
protons along the 23ipa(d, 3n)230u data of Ref. [H is 
covered by the ^^^Pa(p, 2n)^^°U excitation function that 
can be used for the calculation of the inelastic BU en- 
hancement of the (d, 3n) reaction cross sections. On the 
other hand, the following highlighting of the related BU 
effects, discussed at large elsewhere [J-Q, may show that 
there could be different ways to describe the same data 
[l[. While the deuteron-induced fission has previously 
been considered to be the dominant decay channel, the 
same role is attributed in this work to the deuteron BU. 
Consequently, we point out the need of additional mea- 
surements in order to establish which description is bet- 
ter. Nevertheless, a consistent analysis of the deuteron 
interactions is of real interest for applied objectives as the 
fusion technology or nuclear medicine [l[ as well as for 
basic issues related to, e.g., the surrogate nuclear reac- 
tion method (@ and Refs. therein). Thus one may note 
the rising use of the (d, xf) reaction, where x stands for 
a proton or deuteron, as a surrogate for the (n, /) reac- 
tion, and of {d,pj) as a surrogate for neutron capture. 
However, only the internal surrogate ratio method [9| was 



II. DEUTERON BREAKUP AND INDUCED 
FISSION 

A. Deuteron breakup cross sections 

The physical picture of the deuteron breakup in the 
Coulomb and nuclear fields of the target nucleus con- 
siders two distinct processes, namely the elastic breakup 
(EB) in which the target nucleus remains in its ground 
state and none of the deuteron constituents interacts with 
it, and the inelastic breakup or breakup fusion (BF), 
where one of these deuteron constituents interacts with 
the target nucleus while the remaining one is detected 
(e.g., Q and Refs. therein). Under the assumption that 
the inelastic-breakup cross section for neutron emission 
a^^p is the same as that for the proton emission a^p 
(e.g., Ref. [13]), the total breakup cross sections cr_B[/is 
given by the sum 



i/p 

BF 



(1) 



while the total neutron- and proton-emission breakup 
cross sections, a^jj and respectively cr^jj, are given by 



n/p 
C^BU - '^EB 



n/p 
'BF 



(2) 



On the other hand, empirical parameterizations have 
been established [3] for the total nucleon-emission 
breakup fraction 



fin/p) _ n/p I 

Jbu -^buI^R 
and the elastic-breakup fraction 

Ieb = (7eb/(7r. 



(3) 



(4) 
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where api is the deuteron total reaction cross section. 
Since a dependence of these fractions on atomic Z and 
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mass A numbers of the target nucleus, and deuteron inci- 
dent energy E was found on the basis of experimental sys- 
tematics [3|, the nucleon inelastic-breakup fraction may 
have the following form; 



r{n/p) _ r{n/p) r 

JbF - JbU ~ JEB 



(5) 



that leads to the nucleon inelastic-breakup cross sections: 



n/p _ An/p) 



'BF — JBF 



Jbf ■ 



(6) 



However, for ^'^230 only the parametrization of the 
total nucleon-emission breakup fraction [3] 

/^"j/p) = 0.087 - 0.0066Z + O.OOIGSZA^/^ + 

0.0017^^/^£; - 0.000002Z£;2, (7) 

could be accurate, due to the lack of related data for the 
elastic-breakup. In order to use the same approach for 
deuterons incident on the ^'^^Pa target nucleus and the 
energy range of the data of Ref . [l| , we have taken into 
account the value ctbs^IIO mb that was considered for 
deuterons incident on the ^'^^Th nucleus at £^=15 MeV 
[m as well as its normalization corresponding to the mea- 
sured proton inelastic-breakup cross section within the 
same work. Moreover, due to the lack of additional sim- 
ilar data, we used the resulted value /bs =0.164 in the 
whole incident energy range, that is anyway rather nar- 
row. Then we used the corresponding BU and BF com- 
ponents, given by Eqs. ([7]) and (O, to obtain the nucleon 
BF as well as the total BU cross sections shown in Fig. [T] 
In this respect a particular importance has the deuteron 



total reaction cross section in Eqs. ^ and ([6]), respec- 
tively, taken in the present work according to the RIPL-3 
recommendation [l2| for the deuteron optical potential of 
Ref. [l^. Actually this potential is the only one based 
on the data analysis for nuclei with A>208. The corre- 
sponding larger total reaction cross sections (Fig.[T]) with 
respect to, e.g., the default deuteron potential within the 
computer code TALYS led to BU cross sections in 
the present work that are higher than the predictions of 
the Kalbach parameterization [l5| for the total neutron- 
and proton-emission breakup cross-sections that differ by 
normalization factors K, 



d^p/vi 



n/p 
'BU 



K. 



d,n/p 



{A^'^ + 0.8)^ 

1 , (13--E) 

1 -1- exp- — g — '- 



18, K, 



d,p 



21 



(8) 

that are also shown in Fig.[TJ It results that, for deuteron 
incident energies above ~13 MeV, the predictions for the 
total nucleon-emission breakup cross sections given by 
both parameterizations d, [ll| are around ^^50% of the 
deuteron total reaction cross section. It is only the ex- 
trapolation of the Kalbach parameterization at low in- 
cident energies that leads to nucleon-emission BU cross 
sections exceeding even the deuteron total reaction cross 
section. Nevertheless, regardless of the differences be- 
tween them, both parameterizations point out the domi- 
nance of the breakup mechanism at the deuteron incident 
energies below and around the Culomb barrier. Actu- 
ally this conclusion is in line with the experimental total 
proton-emission BU fraction data for deuterons on ^'^^Th 
[111 . [TgI that were taken into account within the above- 
mentioned systematics [3]. 
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FIG. 1: (Color online) Energy dependence of the to- 
tal breakup (solid curve), total nucleon-emission breakup 
cross sections (dashed) [3], and total proton- (dash-dotted 
curve) and neutron-emission breakup cross sections (dash- 
dot-dotted) of Ref. [H, for deuterons interacting with ^^^Pa 
around the Coulomb barrier. The deuteron total reaction 
cross sections [l3l ] are also shown (dotted curve). 



B. Fission competitive decay 

On the other hand, Morgenstern et al. [J] have found 
a definite dominance of the fission decay channel within 
their former analysis of the ^'^^Pa(c?, 37i)^'^°U reaction 
cross-section around the Coulomb barrier. They also 
noted, although without a quantitative assessment, that 
a significant decrease of the available compound-nucleus 
cross section occurs due to the deuteron breakup. Nev- 
ertheless, the fission cross section obtained within the 
EMPIRE-2 computer code assumptions [l3| has been 
quite close to the deuteron total reaction cross section 
(Fig. 3 of Ref. [l|). Conversely, lower fission cross sec- 
tions can be found straight away either using the code 
TALYS- 1.2 with its defauh options [Ti] or within the 
TENDL-2011 library [B], explicitly shown in Fig. [21 The 
two distinct sets of calculated results are used here just 
to prove the similar weight of the fission mechanism. The 
significant difference between them is due to the differ- 
ent optical potentials used for deuterons. Unfortunately, 
there are no (d, /) measured cross sections but only an 
analysis of ^■^^Pa(d,p/) data at the deuteron energy of 15 
MeV p^ . where the measured fission probability versus 
the excitation energy of the residual nucleus ^'^^Pa was 
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FIG. 2: (Color online) Comparison of the excitation func- 
tion of the reaction ^'^^Pa(d,3n)^^''U measured calcu- 
lated using the default options of the code TALYS (solid 
curve), and most recently evaluated [31 (dashed). There are 
also shown the corresponding calculated (dash-dotted) and 
evaluated (dash-dot-dotted) cross sections for the deuteron- 
induced fission on ^^^Pa, together with the evaluated deuteron 
total reaction cross sections (dotted). 



found lower than 40% (Fig. 16 from Ref. 19]). As a 
consequence, we focus on the measured ^^^Pa(d,3n)^"^°U 
excitation function in order to check the dominance of 
the breakup mechanism predicted by empirical parame- 
terizations [sl. [Tsj. 



III. DETAILED BREAKUP MECHANISM 

The dominance of the breakup mechanism has two op- 
posite effects. Firstly, the deuteron total breakup cross 
section reduces significantly the amount of total reaction 
cross section that should be shared among different out- 
going channels. This effect is shown in Fig. EJb) for the 
^'^^Pa(d, 3n)^'^"U reaction and the calculated cross sec- 
tions using the computer code TALYS-1.2 [ij]. In or- 
der to emphasize the two distinct BU effects, we made 
a different choice than for the results shown in Fig. [51 
Subsequently we have not used the BU-inclusion op- 
tion of TALYS by means of the Kalbach parameteri- 
zation [l^. Thus, we have obtained firstly the pre- 
equilibrium (PE) and compound nucleus (CN) contribu- 
tions to the (d, 3n) reaction cross sections, under the as- 
sumption of no breakup process. Then the BU reduction 
of these results was addressed by using a reduction factor 
(1 — <^buI<^r) of the deuteron total reaction cross sec- 
tion. The (d, "in) reaction cross sections obtained in this 
way are now in good agreement with the measured data 
just above the effective reaction threshold while formerly 
these data were also greatly overestimated. However, an 
underestimation by a factor up to 3 at E^IQ MeV be- 



comes visible in Fig. |3l^b). Nevertheless, the description 
of the reaction cross sections at the lowest energies may 
validate the PE and CN model parameters used in these 
calculations. Apart from the default TALYS parameters, 
we used the above-mentioned deuteron optical pot ential 
[Tsl , the nucleon optical potentials for actinides [20[| fr om 
the RIPL 2408 and RIPL 5408 potential segments[il| for 
neutrons and protons, respectively, the microscopic level 
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FIG. 3: (Color online) (a) Results (solid curves) of the con- 
volution of the cross section ratio cr{p,2n) / '^{p,R) for the tar- 
get nucleus ^"^^Pa (dashed) and the Gaussian distribution 
(dotted) of breakup-protons energies for deuterons on ^^^Pa 
at incident energies of 10, 15 and 20 MeV noted on their 
top; in the insert: centroid of the Gaussian distribution of 
breakup-protons energies [2^ versus the deuteron incident 
energy (solid curve) on ^'^^Pa, and the related _Ep ± r/2 val- 
ues (dashed), (b) As for Fig. [5] but for deuteron total re- 
action cross sections [l^l (dotted), nucleon inelastic-breakup 
cross section (short dotted), inelastic BU enhancement (dash- 
dotted), and the (d, 3n) reaction cross sections calculated 
without (dash-dot-dotted) and with (dashed) BU decrease of 
(jR, as well as including the BU enhancement (solid). 
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densities of Hilaire et al. [2l| , and the WKB approxima- 
tion for the fission path model [22j . 

Secondly, we aim to account for the inelastic breakup 
enhancement due to one of the deuteron constituents that 
interacts with the target nucleus and leads to a secondary 
composite nucleus, with further significant contributions 
to various primary reaction channels. The secondary- 
and third-chance particle emissions next to the original 
deuteron-target interaction d, 0, 0| are thus especially 
enhanced. In the present case, the absorbed proton fol- 
lowing the breakup neutron emission contributes to the 
enhancement of the ^aoy activation cross section through 
the ^^^Pa(p, 2n)^^°U reaction. 

In order to calculate the breakup enhancement for the 
^■^^Pa(d, 3rt)^'^''U reaction, the nucleon inelastic-breakup 
cross section iy%p given by Eq. ([B]) was formerly con- 
sidered together with the ratio cf(p^2n)/(^(p.R) that 
corresponds to the weight of the above-mentioned re- 
action induced by the breakup-protons on the ^'^^Pa tar- 
get nucleus We used in this respect the measured 
^'^^Pa(p, 2n)^^'^U reaction cross sections and the pro- 
ton total reaction cross section (T(^p.ii) generated by the 
above-mentioned optical potential. We may express this 
ratio as a function of the deuteron incident energy, using 
the recent Kalbach [2^ formula for the center-of-mass 
system centroid of the Gaussian distribution of breakup- 
proton energies. Since the breakup-proton energy range 
of ^9-14 MeV corresponds to the incident energies of 
11-20 MeV of the measured (d, 3ri) excitation function 
0, the <7(p,2n) cross sections values are provided just by 
the measurements, no other reaction model calculations 
being involved. This is why the simultaneous analysis of 
(d, 3n) Q and (p, 2n) experimental excitation functions 
[2| is so useful for the study of the inelastic breakup and 
complementary reaction mechanisms considered for the 
deuteron interactions with nuclei. 

However, a better estimation for the inelastic breakup 
enhancement is given by the convolution of the ratio 
o'(p,2n)/o'(p,fl) and the Gaussian distribution of the 
breakup-proton energies corresponding to a given 
incident deuteron energy [2^. The former as well as the 
latter quantities for three deuteron incident energies are 
shown in Fig. JSTa). The same method has previously 
been applied [6|, |7[ using a former Gaussian line shape 
[isj . The areas of the related convolution results 
correspond to the inelastic-breakup enhancement of the 
(d, 3n) reaction cross sections at the given deuteron en- 
ergies. The energy dependence of this inelastic-breakup 
enhancement of the ^■^^Pa(c?, 3n)^'^°U activation cross 
section is shown in Fig. [3l^b), while the corresponding 



total activation of ^^o-jj jg g^g^yy compared with the 
experimental data As expected, the more realistic 
treatment of the inelastic breakup enhancement by 
taking into account the quite large widths F of the 
breakup-proton energy distributions (see the upper 
insertion in Fig. 3(a)), has led to a rather accurate 
description of data . Further improvements of the 
breakup analysis could add to a better account of the 
related energy dependence. 



IV. CONCLUSIONS 

An analysis of the ^'^^Pa((i, 37i)^'^'^U reaction excita- 
tion function at energies around the Coulomb barrier has 
taken into account the pre-equilibrium and compound- 
nucleus cross sections corrected for the deuteron-breakup 
decrease of the total reaction cross section, as well as 
the inelastic breakup enhancement. The analysis reveals 
the dominance of the deuteron breakup mechanism un- 
like a former assessment in this respect of the deuteron- 
induced fission process. 

The improvement of deuteron breakup effects esti- 
mation requires complementary experimental studies of, 
e.g., the (c?, 3n) and (p, 2n) reaction cross sections for 
the same target nucleus and within correlated incident- 
energy ranges. The suitability of the empirical parame- 
terization of the breakup components and reaction mech- 
anisms involved in the interaction process could thus be 
checked and updated. Furthermore, the associated in- 
clusive neutrons and protons spectra measurements that 
allow the distinction among various contributing mecha- 
nisms are highly requested too, as well as {d,pf) angular 
correlations when the deuteron induced fission process 
is analyzed. Given the increased interest in surrogate 
reaction studies, e.g. Ref. [13], the usefulness of de- 
tailed theoretical and experimental investigations of the 
breakup of weakly bound projectiles including deuterons 
is obvious. 
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